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Jr^ ' In this article we describe recent progress in the computational many-body theory of metal sur- 

Q\ faces, and focus on current techniques beyond the local-density approximation of density-functional 

theory. We overview various applications to ground and excited states. We discuss the exchange- 
' correlation hole, the surface energy, and the work function of jellium surfaces, as obtained within the 

random-phase approximation, a time-dependent density-functional approach, and quantum Monte 
Carlo methods. We also present a survey of recent quasiparticle calculations of unoccupied states 
' at both jellium and real surfaces. 
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I. INTRODUCTION 



Prerequisite for calculating the electronic properties of solid surfaces is a practical procedure for the treatment of 
, electron-electron (e-e) interactions in large, strongly inhomogeneous electron systems. Over the last three decades, 
^ • almost|-|all first-principles electronic-structure calculations have been carried out within density-functional theory 
O ' (DFT)lil, which solves the many-electron problem by introducing an auxiliary non-interacting system with the same 
electron densitycl. Although DFT is a formally exact theory, it relies on the knowledge of the so-called exchange- 
correlation (xc) energy as a functional of the electron density. The best known approximation for this functional is the 
local-density approximatioB, (LDA)0, which has proved remarkably successful in various applications to ground-state 
^ ' properties of solid surfacescl. However, this simple approximation presents two important shortcomings. First, the 
LDA is unable to provide the highly accurate total energies that are required in the application of DFT to chemical 
^ ■ reactions at surfaces. Secondly, it yields an inaccurate description of the long-range xc-hole density at solid surfaces, 
] and fails to reproduce the correct image- like asymptotic behaviour of the surface barrier, which plays an important 
role in the interpretation of a number of surface-sensitive experiments. Furthermore, the knowledge of the finite 
' quasiparticle lifetimes of unoccupied states is beyond the scope of DFT. 

, The aim of this article is to survey recent computational work on the description of many-body effects in both 
ground- and excited-state properties of metal surfaces. In section II we give a short overview of the state-of-the-art 
computational many-body theories of the inhomogeneous electron gas: DFT, in the LDA and its various gradient- 
corrected forms, quantum Monte Carlo (QMC) methods, and Green function theory in the so-called GW approxima- 
tion. In section III we focus on a description of ground-state properties of metal surfaces, and we critically examine 
the performance of current many-body schemes in the computation of two of the most fundamental magnitudes, 
Q I namely the surface energy and the work function. In section IV we examine quasiparticle calculations of the, surface 
O barrier and the binding energies and lifetimes of image-potential induced states. These so-called image stateaJ form a 
Rydberg-like series, which converges towards the vacuum energy, and their existence is perhaps the most striking man- 
ifestation of long-range many-body effects at surfaces. We conclude this article in section V with some perspectives 
for future research in this area. 

We use Hartree atomic units throughout, unless otherwise specified. In these units, e — rUe — fi — 4,Treo — 1. 



II. THEORY OF THE INHOMOGENEOUS ELECTRON GAS 

Under the assumption that electronic and ionic degrees of freedom can be decoupled, the problem of N electrons 
in a solid can be described by the following Hamiltonian: 

" r 1 1^1 

^ = E -2^? + K..(rO +Et^-H- (2-1) 

1=1 L J i<j ' * 
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Here is the coordinate of electron i, and Vextif^i) is the external potential accounting for the interaction with the 
fixed atomic nuclei. 

A. Density-functional theory 

DFT is a formally exact theory based on two theorems by Hohenberg and Kohrfl. The first theorem is extraordinary: 
it states that the ground-state density n{r) of a system of interacting electrons determines the ground-state properties 
of the system uniquely. The second theorem states that the ground-state energy can be obtained by minimising an 
energy functional i?[n(r)] of the electron density. The minimum value of i?[r;,(r)] is the exact ground-state energy, 
and is attained when h{r) is the exact gMund-state density n(r). 

In the Kohn-Sham formulation of DFTu, the ground-state density can be obtained by calculating the eigenfunctions 
4>i(r) and eigenvalues Si of non-interacting, single-particle equations, yielding 



n(r) 



N 

E 

i=l 



(2.2) 



and the energy functional is written in the form 



E[n{r)] = To[n(r)] + / drn(r)K.t(r) + - 



, n(r)n(r') ^ r~/ ^^ 



(2.3) 



Here To[n(r)] is the kinetic-energy functional for non-interacting electrons with density n{r), and Exc[n{r)] is the 
so-called xc energy functional. E^c contains all many-body effects beyond the Hartree approximation, as well as 
the difference between the ac tual kinetic-energy functional of the interacting system and To[n(r)]. By minimising 
the energy functional of Eq. ( 2^ ) , one concludes that the exact density is obtained by solving self-consistently the 
so-called Kohn-Sham equations 



-y + Vextir) + 



(2.4) 



where 14c (r) is the xc potential. 



vu^) 



5Exc[n{r)] 



(5n(r) 



(2.5) 



n(r)— n(r) 



The xc energy functional is not known and must be approximated. The simplest approximation for Exc[n{'ry\ is the 
LDA, 



£;L°^[n(r)]= j drn(r)C'-^(n(r)), 



(2.6) 



where the actual xc energy per particle at point r, ea;c(r; [n(f)]), which is a functional of the density n(f), has been 
replaced by the xc energy per particle of a uniform electron gas of density n(r), e""*-^(n(r))[]. 

For an electron gas of slowly varying density the LDA may be corrected through the introduction of gradient 
expansions, as suggested by Kohn and ShamEl. Although the basic physics behind gradient expansions is correct, 
Langreth and Perdewl3 showed that it fails in most practical situations because of the spurious small-wave-vector 
contribution to the Fourier transform of the second-order density-gradient expansion for the xc hole around jaii 
electron. Via cutoff of this spurious contribution, generalised gradient approximations (GGA) were then devisedEl^llj, 
which approximate the xc energy per particle as a carefully chosen nonlinear function of the election density and 
its gradient. Recently, a meta-generalised gradient approximation (meta-GGA) has been developedtj, which makes 
use not only of the local density and its gradient but also of the orbital kinetic-energy density. For other kinds of 
improwjiHents of the LDA, including the weighted density approxjimation (WDA)liJ and self-interaction corrections 
(SIC)llall3, we refer the reader to the review by Perdew and Kurthll3. 



^Several parameterisations of e^dn) exist, which are based on many-body QMC simulations of the uniform electron gasi 
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B. Quantum Monte Carlo methods 



Unlike the density-functional approach, in which the ground-state density is the basic quantity, wave-function- 
based quantuta Monte Carlo methods sample the full ground-state many-body wave function of the system under 
considerationES. There are many different QMC methods, but here we will be concerned with only two of them: 
variational quantum Monte Carlo (VMC) and fixed- node diffusion quantum Monte Carlo (DMC), which have been 
shown to provide a description of increasing sophistication of the properties of zero-temperature systems. 

The VMC method is based on an ansatz for a parameterised trial many-body wave function. Expectation values 
are evaluated by the Metropolis Monte Carlo techniquesEj'EEl, and the free parameters are then varied in order to 
optimise either the energy expectation value or the fluctuatiop^of the local energyE3. 

We consider trial wave functions of the Slater- Jastrow typeO 



vl/ = D exp 



N 

E 

1=1 



N 



(2.7) 



where Z? is a Slater determinant of either Hartree-Fock (HF) or LDA single-particle orbital^, Ua--,a-j is a two-body 
term correlating the motion of pairs of electrons, asd Oi denotes the spin of electron i. The short-range form of 
Uai,<7j is dictated by the weU^mpwn cusp conditionE3, while the long-range behaviour is related to the zero-point 
motion of plasma oscillationaH^'Ej. The one-body term Xo-. (r), which is absent in the case of a uniform electron gas, 
allows a variational adjustment of the electron density in the presence of the two-body term. Despite the apparent 
simplicity of the Slater- Jastrow trial wave function of Eq. (|]^), VMC has proved very successful in applicatiq 
bulk properties of solids. In particular, VMC methods were used to calculate cohesive energies of solid materials 
showing that they are typically an order of magnitude more accurate than cohesive energies obtained within the HF 
or LDA approaches. 

The main drawback of VMC is that the accuracy of the results is entirely dependent on the quality of the trial wave 
function. The DMC method overcomes thisJimitation, by using a projection technique to enhance the ground-state 
component of a starting trial wave functionEj. In DMC, the mathematical equivalence between the imaginary-time 
many-electron Schrodinger equation and a real-time diffusion equation is used to solve the former by a simulated 
numerical diffusion of a collection of fictitious particles (walkers). Each walker corresponds to a point in the 3iV 
dimensional space of electron configurations. The density of walkers represents the many-electron wave function, and 
the propagation of walkers, according to the imaginary-time Schrodinger equation, exponentially projects out the 
ground state from the initial state (provided that the initial state is not orthogonal to the ground state) . In this way, 
an exact numerical representation of the ground state can be obtained from a reasonably accurate initial guess. The 
initial wave function for the more accurate and computationally expensive DMC calculation is usually provided by 
VMC methods. 

The analogy between the imaginary-time Schrodinger equation and the diffusion equation only works if the many- 
body wave function is a positive (or negative) function, which can then be interpreted as a density distribution. As 
for fermions the antisymmetric many-electron wave function must take both positive and negative values, there is 
a sign problem in DMC caJpilations for these particles. The common way to circumvent this problem is to use the 
fixed-node approximationEju. Within this approach, for a given trial many-electron wave function one defines a trial 
many-electron nodal surface on which the trial function is zero and across which it changes sign. The fixed-node 
DMC algorithm then produces the lowest-energy many-electron state with the given nodal surface. The resulting 
fixed-node errors are typically about 5% of the correlation energy (defined as the difference between the exact and the 
HF ground-state energies )e2I. In most systems exchange dominates the correlation by almost an order of magnitude, 
and fixed-node errors are, therefore, reduced to about 0.5% of the total ground-state energy. Another possible source 
for systematic error in QMC calculations is the use of a finite simulation cell to model an extended system]^ We 
shall return to these finite-size errors in section III.B, where we will examine QMC calculations of the jellium surface 
energy. 

Although VMC and DMC methods are best suited to study the ground state, they have been shown tOjalso provide 
some information about specific excited states. This subject has recently been reviewed by Foulkes et aZEil. 



is often split into electron spin-up and spin-down Slater determinants, in order to reduce computational cost. 
^Within a many-electron formulation, the Schrodinger equation cannot be solved by considering only a single primitive unit 
cell of the periodic nuclear potential, since the e-e interaction breaks the lattice symmetry. 
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C. Green-function theory 



The Kohn-Sham eigenvalues are often found to be in good agreement with the quasiparticle energies measured in 
photoemission experiments. Nevertheless, they have no clear physical meaning, except for the highest occupied one 
which corresponds to the negative of the ionisation energyH, and there are some serious discrepancies with measured 
single-particle excitation energies. Also, the knowledge of the finite lifetime of excited states is known to be outside 
the domain of DFT. 

A ppaper way of calculating single-particle excitation energies and lifetimes is provided by the Green-function 
theoryEj and the quasiparticle concept accounting, in a many-electron system, for the ensemble of each electron and 
its own xc hole. With |A^, s > as some many-body eigenstate s {s — OJar the ground state) of N electrons, one defines 
the quasiparticle amplitude ^'^(r) and the quasiparticle energy Eg aa^il 



*,(r) =< N,0\,p{r)\N + l,s>, 



E 



N.O 



for Eg > fj. 



(2.8) 



and 



*,(r) =< N -l,.s\^j{r)\N,0>, 



Eg — En^o — E 



for Eg < fi, 



(2.9) 



where ip{v) is a field operator in the Heisenberg representation that destroys an electron at point r, and /i represents 
the chemical potential. |— . 

The one-particle Green functioned determines the energy and the damping of the quasiparticles. This can be shown 
explicitly by introducing the complete set of many-body eigenstates of the full Hamiltonian for the -I- 1 or iV — 1 
particles in the definition of the one-particle Green function. One obtains 



G(r,r';0=-*E*^W*^We~ 



iE.t 



(2.10) 



which identifies —2lniEg as the inverse quasiparticle lifetime. Furthermore, Fourier transformation of Eq. (2.10) to 
energy space yields 



Gir,r';E) = J2 

S 

which is known to satisfy the Dyson equation 



^s{r)^*g{r') 
E-Es ' 



G{r,r';E)^G"{r,r';E)+ J dr, J dr^ G°(r, n; i?) S(ri, ra; i?) G(r2, r'; i?). 



(2.11) 



(2.12) 



Here G"(r, r'; E) represents the non-interacting one-particle Green function, and S(r, r'; E) is the so-called self-energy 
of the quasiparticle, which is a nonlocal, energy-dependent, non-Hermitian operator accounting for all xc effects 
beyond the Hartree approximation. By inserting the spectral representation of Eq. (2.11) into Eq. ( 2.12] ), one finds 
the quasiparticle equation 



ly' + Ve.tir)+ J dr'j^^ ^g{r)+ J dr' J:{r,r'; Eg)^gir') = Eg^g{r). 



(2.13) 



For a system of interacting electrons, there is little hope in solving the quasiparticle equation exactly. Hence, one 
usually resorts to perturbation theory and chooses, as a starting point, a suitable single-particle Hamiltonian Hq 
whose wave functions </'s(r) and energies Eg approach the quasiparticle wave functions ^'s(r) and energies Eg. In usual 
practice, one takes the LDA Kohn-Sham Hamiltonian 



LDA 



(r), 



(2.14) 



where 



LDA 



(r) 



d[nejf™-^(n)] 



dn 



(2.15) 



n— n(r) 
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The quasiparticle energy is then determined from 

Es-e 

where 



and 



AE,{E)^ Jdr Jdr'^sir) [^{r,r' ; E) ~ VUr)6{r - r')] M^'), 

dAJ:s{E) 



dE 



E=E, 



(2.16) 
(2.17) 

(2.18) 



is the so-called quasiparticle weight or renormalisation factor. On the energy-shell, this factor is takep to be unity. 

The exact self-energy can be obtained, in principle, from an iterative solution of Hedin's equationa^J in combination 
with the Dyson equation. However, to obtain explicit results one usually resorts to an expansion in powers of the time- 
ordered screened interaction W{r, r';E). The leading-order term of this expansion is the so-called GW approximation: 



— e-'">^' G(r, r'; E - E') W{r, r'; E'), 
-oo 27r 



(2.19) 



which can also be obtained as the first iteration of Hedin's equations by simply neglecting vertex corrections. The 
screened interaction can be expressed in terms of the density-response function x(r, r'; E), as follows 



W{r,r';E) ^v{r-r') + J dn J dra w(r - n) x(ri, ra, S) w(r2 - r'), 



(2.20) 



where w(r — r') repr esents the bare Coulomb interaction. If one replaces the screened interaction by w(r — r') the self- 
energy of Eq. ( ^.19 ) reduces to the Hartree-Fock self-energy, which sometimes leads to unphysical results, especially 

in the case of metals where screening plays a crucial role. 

Most current GW calculations simply replace the exact one-particle Green function entering Eq. (|2.19| ) by the 
non-interacting Green function G°{r,r'; E) pertaining to the LDA Kohn-Sham Hamiltonian of Eq. (^TJ)7 On the 
same level of approximation and neglecting all vertex corrections, the screened interaction entering Eq. (2.19) is 
usually obtained from Eq. ( ^.20| ) with the density-response function evaluated in the random-phase approximation 
(RPA)g: 

X(r,r';i?)=x°(r,r';i?) + J dn J dra x°(r, n; i?) «(ri - ra) x(r2, r'; £;), (2.21) 

X°(r, r';^^) being the non-interacting density-response function 

X"(r, r'; E) ^ J dE' G"(r, r'; E') G°{r, r'; E + E'). (2.22) 

A survey of the theory underlying the GW formalism and its applications can be found in Refs. |3^ and |3^. 
For improvements of the standard [G^W^) quasiparticle calculations, including self-consistency and beyond-GW 
approximations, we refer the reader to the recent review by Aulbur, Jonsson, and Wilking22l. 

III. GROUND-STATE PROPERTIES OF SURFACES 

The simplest model to investigate the many-body properties of simple metal surfaces is the well-known jellium 
model. Within this model, valence electrons are described by an inhomogeneous assembly of free electrons moving 
in a uniform neutralising positive background, which is abruptly truncated at a plane. In this section, we focus on a 
description of many-body properties of jellium surfaces, and we also survey, when available, more realistic calculations 
where the effect of the discrete crystal lattice is taken into account. In the discussion of the jellium surface, the z axis 
is taken to be normal to the surface, with the metal occupying the z > half-space. The positive-background charge 
density n is expressed in terms of the Wigner radius r^, as 1/n = (47r/3)r^. 



This DFT-based RPA differs from the less realist ic ac tual RPA, which is determined from the non-interacting Green function 



associated with the Hartree [V^c (r) = in Eq. (2.14)] Hamiltonian 
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A. Pair-distribution function and exchange-correlation hole 



Central quantities in the discussion of xc effects in an interacting many-electron system are the pair-distribution 
function ^(r, r') and the xc-hole charge density n2;c(r,r'). The pair-distribution function is defined as the probability 
of finding an electron at point r' if there is already one electron at r 



<?(r,r') = 



n(r)n(r' 



- [< N,0\h{r)n{r')\N,0 > -n{r)S{r - r')] 



(3.1) 



where \N, > represents the many-body ground state of N interacting electrons, n(r) is the electron density, apd "•(!■) 
represents the electron-density operator. Making use of the zero-temperature fluctuation-dissipation theorerrH, and 
by a suitable choice of integration contoursEj, the pair correlation function can be written in terms of the retarded 
density-response function x(r, r'; E) of the many-electron system^ 



5(r,r') = l 



1 



1 r°° 

- dEx{r,r';iE)-n{r)S{r-r') 
Jo 



n(r)n(r') 

The total interaction energy Eee of the many-electron system is obtained as follows 

2. I r-r' 



which can be written in the form 



Eee = J drn{r)eeeir; [n{r)]), 



(3.2) 



(3.3) 



(3.4) 



eee(r; ['T-(f)]) representing the interaction energy per particle at point r. This interaction energy can be separated into 
electrostatic (Hartree) and xc terms, 



eee(r; [n(f)]) 



dr' 



, n(r') 



1 



dr' 



, n^c(r, r') 



|r-r'| 2, 

na;c(i", r') being the so-called xc-hole charge density 

n,,(r,r')-n(r')[l-5(r,r')] 



r — r' 



(3.5) 



(3.6) 



This quantity represents the depletion of the average electron density n(r') due to the presence of an electron at r, 
which is dictated by the combined effect of the Pauli exclusion principle and Coulomb correlations, and satisfies the 
sum- rule 



J dr' n^e{r,r') = 1. 



(3.7) 



Eqs. ( pTl| ) pi|d-|(3.6) represented the starting point for raeent QMC calculation s o f n^c in bulk siliconH, the 
jellium surfacecilEHl, and strongly inhomogeneous model solidald. Alternatively, Eq. (3^) allows a calculation of the 
pair-distribution function and the xc hole from the knowledg e of the density-response function of the system. In a 
DFT-based RPA, this function is usually obtained from Eq. (2.21) with the use of the single-particle eigenfunctions 
and eigenvalues (j^fiAc LDA Kohn-Sham Hamiltonian of Eq. (2.14). In the framework of time-dependent density- 
functional theoryESca, the exact density-response function satisfies the integral equation 



X(r,r';£;) =x°(r,r';i;) + J dn J dra x"ir,ri;E) [«(ri, ra) /.c(ri, rs; i;)] x(r2, r'; £;), 



(3.8) 



where 



^For positive energies this function coincides with its time-ordered counterpart, which enters the self-energy formalism de- 
scribed in the preceding section. 
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/..(r,r^i^) = %|!^, (3.9) 

Vxc [n{r, E)] representing the exact energy-dependent xc potential, and x^i^i ^'i being the density-response function 
of non- interacting Kohn-Sham electrons. In the so-called time-dependent local-density approximation (TDLDA)La 
or, equivalently, adiabatic local-densi ty a pproximation (ALDA), the Kohn-Sham effective potential is obtained in the 
LDA and the exact xc kernel of Eq. (|3.9|) is replaced, for all energies, by 



5{r-v'). (3.10) 

n— n(r) 



Assuming, within a jellium model of the surface, translational invariance in the plane normal to the z axis, the 
xc-hole charge density is found to be given by the following expression: 



nxc{r;z,z') = I dq^^ Jo{q\\r) 

2Trn(z) ' 



1 



dE xiz, z'; qii ,iE)]+ n{z)6{z - z) 



(3.11) 



which depends on the distance r — \vn — ry| projected in the plane of the surface. If the actual density-response 



function x{'^t'^'\%t^^^ is replaced by its non-interacting counterpart ^(z^z'\q\^,iE\ then Eq. (3.11) yields the 
exchange-only or Fermi hole^. 

In the simplest possible microscopic model of the jellium surface, namely the so-called infinite barrier model (IBM), 
the one-electron wave functions describing motion nornaal to the surface are simply sine functions. Using these wave 
functions, the exchange hole was studied by JuretschkecZl and otherscS, and the xc hole was later investigated in the 
RPA by Inglesfield and MooreEJ. However, the IBM does not permit the electron density to relax beyond the infinitely 
high potential barrier outside the surface, and yields a poor representation of both the electron density and the xc 
hole. 

The behavipur of the exchange hole for an electron that is located well outside the surface was investigated by Sahni 
and BohnenEj with the use of one-electron wave functions that were generated from a linear-potential model. These 
authors found that inside the solid, and up to the position of the jellium edge, the x-hole density rixij"., z, z') shows a 
behaviour similar to that obtained within the IBM. As the electron overcomes the jellium edge, the exchange hole was 
found to remain behind and to be localised near the surface region, as within the IBM. However, when the electron is 
removed further away from the surface into the vacuum the exchange hole, instead of staying localised at the surface, 
was found to spread ttmoughout the entire solid. This behaviour was further confirmed by the calculations reported 
by Harbola and SahniEil, who used a step-potential model of the jellium surface. These authors also considered the 
x-hole charge distribution in the plane parallel to the surface, and found that it differs significantly from the classical 
image charge distribution. 

Recently, preliminary self-consistent calculaticps of both exchange and corr elatio n co ntributions to the xc-hole 



charge density at a jellium surface were reportedEj, as obtained from Eqs. (3_^) and (3^) with a full description of 
the DFT-based RPA density-response function of a jellium slab. These calculations showed that for an electron well 
outside the surface the x-only charge density displays several maxima within the solid, spaced at intervals of half a 
Fermi wavelength, in agreement with the results of Sahni and Bohnen, and confirmed the interesting, albeit expected 
result, that the build up of the image charge is entirely due to the surface-localised part of the correlation hole. 

Fig. 1 shows contour plots of the fully self-consistent DFT-based RPA xc hole near a jellium surface with the 
positive-background charge density n of aluminium (r^ = 2.07). In the interior of the metal the xc hole is a sphere 
cantered at the position of the electron, as in the case of a uniform electron gas. As the electron moves towards the 
jellium edge, the xc-hole charge density n^c begins to distort and starts to lag behind. Finally, as the electron is 
removed far outside the surface the xc hole flattens and remains localised near the surface, becoming an image-charge 
distribution located at the so-called image plane. 

The issue of the physical origin of the image potential has been the source of considerable controversy over the 
years. The asymptptic behaviour of the x(j-^otential V^c(r) of DFT at large distances outside a metal surface was 
examined by ShamE3 and by Eguiluz et a/O through the following exact integral equation, which relates Vxci^) to 
the xc self-energy E2;c(r, r'; i?)P 



(z,z' \q\\,iE) were determined from the eigenfunctions and eigenvalues of the non-local Hartree-Fock Hamiltonian, this 
Fermi hole would coincide with that obtained in the Hartree-Fock approximation. 
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J dr'T4e(r') J dEG°{r,r';E)G{r',r;E) = J dr, J dv2 J dE G''{r,r,; E) S,,(ri, ra; £;) G(r2, r'; i;). 



(3.12) 



Both ShanH and Eguiluz et alE3 concluded that the exchange self-energy Sa;(r,r';i?) yields a 1/z^ behaviour for 
large z, while the correlation self-energy yields the image-like 1/z limit, thereby confirming the long-standing believe 
that the image-potential structure is a pure Coulomb-correlation effect. 

Harbola and Salwii suggested an alternative procedure for constructing the xc potential that enters the Kohn-Sham 
equations of DFTE3. These authors interpreted 14c (r) as the work required to remove an electron from the solid 
against the electric field of its xc-holc charge distribution, and showed that the exchange-only potential 14 (r) merges 
with the classical image potential. Although this conclusion seemed to be in disagreement with the numerical jellium- 
slab G^VF^ calculations reported in Ref. it was then argued that while correlation gives the asymptotic image 
limit for slab geometries exchange, does provide the limit in the case of the semi-infinite solidEJ. Nevertheless, the 
conclusion of Harbola and SahniE3 is still in disagreement with the earlier investigations of Ref. |5^, and this issue 
deserves further investigation. Alternatively, the image potential felt by electrons well outside the surface can be 



obtained from the asymptotic behaviour of the interaction energy per particle eee(r; ['T-(f)]) of Eq. (3.5), which, is 



dictated by the interaction between each electron and its own xc hole. Work along these lines is now in progress 
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B. Surface energy 

One of the surface properties most sensitive to many-body effects is the surface energy. This static physical quantity 
is defined as the energy per unit area required to split the solid into two separate halves along a plane, i.e., 

2E[n{r)]-E[n\r)] 
o= , (3.13) 

where n(y) represents the ground-state electron-density distribution for each half of the solid after it is split, and n'(r) 
is that for the unsplit solid. Following the DPT energy-functional partitioning of Eq. (pT^), the surface energy can be 
written as the sum of three terms: 

Cr = (Ts + (Jes + (Ja;c, (3.14) 

where CTs represents the kinetic surface energy of non-interacting Kohn-Sham electrons, Ues is the electrostatic surface 
energy due to all positive and negative charge distributions in the solid, and a^c is the xc contribution to the surface 
energy. For a jellium surface, one writes 

dzn{z) {e,e(z; Hz)]) - e^T^in)} , (3.15) 



where £xciz; [n{z)]) is the xc energy per particle at point z and £"™-^(n) is the xc energy per particle of a uniform 
electron gas of density n. In the LDA, Exdz; [n{z)]) is simply replaced by e""*-^(n(z)). I-. 

The first self-consistent DFT calculation of the jellium surface energy was carried out by Lang and KohncI, in the 
LDA approximation. Lang and Kohn also discussed the effect of the crystal lattice, within first-cmier perturbation 
theory, and this work was later generalised to treat the discrete-lattice perturbation variationalljO. Nowadays, the 
full inclusion, within the LDA, of lattice effects poses no difficulties; however, the question of the impact of nonlocal 
xc effects on the surface energy, and their interplay with the strong charge inhomogeneity at the surface, has remained 
unsettled over the years. 

The original discussion about the effecL-of nonlocal Coulomb corpjations was centred around the surface-plasmon 
contribution to the surface energy. CraigjEl and Schmidt and LucatO independently claimed that the surface energy 
of simple metals can be accounted for by the change of zero-point energy of the plasma oscillations resulting from the 
cleavage of the metal, which is a nonlocal effect not contained in the LDA-calculation of Lang and Kohn. However, 
the actual magnitude of this contribution became a matter of controversyEa, and more careful calculations were then 
presented. i-Ijhe xc surface energy was first obtained in the RPA for a free-electron gas that was terminated abruptly 
at a planeE^'EJ, and more elaboraie RPA calculations were then reported for the case in which the surface potential 
is taken to be an infinite barrierS A comparison between these IBM-based RPA calculations and a local-density 
calculation of <Jxc for the IBM density profile yielded a difference-pf ~ 10%, though much larger differences were found 
for the separate exchange and correlation contributions to (T2:cE3. Lang and Sham then argued that a regular power 
series in density gradients does not exist for separate exchange and correlation energies, apd showed that the LDA is 
satisfactory for the combined xc term but not for exchange and for correlation separatelyEZI. 
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RPA and ALDA calculations of jellium surface energies that are based on a fulUin self-consistent description of 
the one-particle Hartree or LDA orbitals have been performed only very recentlyE3E3. In these calculations, both 
exc{z] [niz)]) and e^J-Ln) were obtained from a coupling-constant integration of the xc contribution to the interaction 

£.e(r; [n(f)]) = \ d\ j dv' (3.16) 



energy of Eq. (p 



where n^^{Y,Y') is the xc-hole charge density of a fictitious system so chosen that as the e-e interaction strength A is 
varied in the interval [0, 1] the corresponding density equals its fully interacting value, i.e., 'n\{v) = n[r). In Refs. ^ 
and ^ LDA surface energies were also computed, as obtained from the RPA and ALDA xc energies per particle of 
a uniform electron gas, and a comparison of local versus nonlocal surface energies supported the conclusion that the 
error introduced by the local-density approximation is small. 

The RPA is exact for exchange and long-range rGcrrelations, but it is known to be a poor approximation for 
short-range correlations. Hence, Kurth and Perdewllj constructed a GGA for the short-range correlation energy in 
combination with the DFT-based RPA calculations of Refs. |6^ and ^ for long-range correlation, and obtained xc 
surface energies that are larger than the LDA prediction by about 2%. Furthermore, it was argued that while the 
RPA is not a good approximation for the total correlation energy, it seems to be a S|Wfprisingly good approximation 
for certain changes in the correlation energy, such as those arising in surface energieaI3. 

Two new approaches have been developed recently, which yield sucface energies agreeing to within 1% with the 
so-called RPA -I- calculations of Ref. These are: (a) a metaJISIGAEj, and (b) a new wave-vector interpolation as 
a long-range correction to the GGA exchange and correlatioro. Nonetheless, thesfi-jCalculations strongly disagree 
with the result of wave- function-based Fermi- hypernetted-chain (FIINC)Ej and DMCcda calculations, which predict 
surface energies that axe significantly higher than those obtained in the local-density approximation. 

Acioli and CeperleyCJ performed fixed-node DMC calculations of the surface energy for a jellium slab that is con- 
tained in a finite simulation cell. They obtained surface energies that are systematically larger than the corresponding 
LDA results. Although for a high-density metal (r^ — 2.07) the DMC surface energy was found to be close to the LDA 
prediction, the difference with the LDA surface eiiergy was found to be as large as 50% for = 4, in agreement with 
the FHNC calculations of Krotscheck and KohnLJ (the absolute difference between DMC surface energies and those 
obtained in the LDA was found to be of about 150 erg/cm^ for all densities under study, i.e., 2.07 < < 4). The 
discrepancy between DFT and wave-function-based surface energies is surprising: while fixed-node DMC | CalC | Yil at4( 
are often regarded as essentially exact, all the current DFT-based estimates of the jellium surface energyliSHoO 
seem to be consistent with the conclusion that LDA-and GGA calculations for real surfaces provide energies at least 
as accurate as those derived from the experimentiljilS. 

A possible source of systematic error in QMC calculations for extended systems is the use of a finite simulation 
cell. Hence, one needs to extrapolate the ground-state energy Eq^ of the extended medium from that obtained for 
the A'^-particle system, Ef^. In the case of a uniform system, one usually uses the extrapolation formula 

E^^EN+a{T^-TN) + ^, (3.17) 

where the parameters a and h are fitted by considering the energy of increasingly larger systems. Too and Tjv represent 
the kinetic energies of the non-interacting infinite and finite systems, and the b/N term accounts for the so-called 
Coulomb finite-size effects arising from the e-e interaction-energy term. The non-interacting term (Too — T/v) is 
generally obtained with the use of one-particle LDA calculations, and the parameter a (a ~ 1) accounts for the 
difference between the kinetic energies of the interacting and the non-interacting systems. For non-uniform systems, 
there are also finite-size errors due to the electrostatic-energy contribution to the ground-state energy, which only 
depend on the electron density. As the LDA densities are found to be rather accurate, these finite-size errors can also 
be estimated within tlie LDA. 

Acioli and Ceperljt2l estimated the finite-size corrections to their DMC surface energies from the difference between 
the LDA surface energies of their finite simulation cell and a jellium slab that is infinitely extended in the plane 
parallel to the surface. Although this estimate accounts for kinetic and (part of) electrostatic contributions to the 
finite-size effects, there are additional (Coulomb) finite-size errors, not encountered in the LDA, which are originated 
in the periodic Ewald interaction used to model the e-e Coulomb interaction in a periodic geometry. 

In order to estimate the magnitude of the Coulomb finite-size corrections, which are not included in the DMC 
surface energies of Ref. ^ we used the VMC method to calculate the surface ene rgy of a jellium slab with = 3.25. 
For the trial wave function we, used the standard Slater-Jastrow form of Eq. ( ^.7[ ) in combination with the one-body 
term Xi (r) due to Fahy et alx3. The calculations were performed with the use of a faced-cantered cubic (fee) supercell 
containing N = 284 spin-unpolarised electrons. We modelled the e-e interaction energy using either the standard 
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Ewald interaction, as used by Acioli and Ceperleyc^, or the model periodic Coulomb (MPC) interaction of Ref. |8^, 
which virtually eliminates Coulomb finite-size effects. Our calculations indicate that for N = 284 and Vg = 3.25 
Coulomb finite-size corrections to the VMC simulations are positive and amount to 3.9% of the surface energy. Since 
these corrections decay as with the size of the system, we estimate that Coulomb finite-size corrections to the 
DMC calculations of Acioli and Ceperley amount to ~ 6% of the jellium surface energy for — 3.25. Our preliminary 
VMC calculations yield surface energies of 433 and 450 erg/cm^, respectively, depending on whether the Ewald or the 
MPC interaction is employed. These surface energies are higher than predicted by the DMC calculations of Acioli 
and Ceperley, which yield for Vg = 3.25 a surface energy of 360 erg/cm^, and are well above the LDA prediction of 
220 erg/cm^ for this electron densityj]. 

Hence, finite-size corrections do not seem to account for the existing large differences between QMC and LDA surface 
energies, which are not present in DFT-based calculations, and a satisfactory understanding of nonlocal correlation 
effects in the surface energy of simple metals remains to be given in the future. Although the LDA and GGA for real 
metals have been found to provide surface energies in fair agreement with the experiment, this might well be due to 
a cancellation between the errors introduced by the LDA in the total energy of both the infinite and the semi-infinite 
solid. Such a cancellation of errors could not occur in the case of jellium surface energies, as the LDA is exact for the 
infinitely extended jellium. 



C. Work function 



The work function $ is the minimum energy required to remove an electron from the solid to a macroscopic distance 
outside the surface. In a many-body context (and at T = 0) this corresponds to a process in which a many-body 
system composed of N electrons is brought from its ground state to an excited state with — 1 electrons in the lowest 
possible state in the metal and one electron at rest at infinity. Hence, one writes 



<^> = En' - En, 



(3.18) 



where En is the ground-state energy of the semi-infinite solid and En' is the energy of the excited state. The long- 
range image potential tail of the surface-electron interaction potential ensures that the electron remains bounded to 
the surface, occupying a bound state n (n = 1 for the lowest-lying state) of the image potential with encrgyQ 



1 



4(n + a)2 



(3.19) 



where a is the so-called quantum defect, and 0oo is the electrostatic potential well outside the surface. Thus, the 
minimum energy of an electron moving freely at a macroscopic distance outside the surface is 



= lim 



1 



4(n -I- ay 



and the total energy of the lowest-lying excited state of the semi-infinite solid with one electron at infinity is 



E' 



N 



En-1 + I 



(3.20) 



(3.21) 



As the highest occupied DFT eigenvalue Ep equals the (negative of) the ionisation energy En — En-\, Eq. ( 3. IS ) 
may be written in the equivalent one-body form 



Ef 



(3.22) 



or, equivalently, 



$ = 1:^4)- Ep, 



(3.23) 



'^In these calculations we used the DMC total energy of as quoted by Acioli and Ceperley for the total energy of the uniform 
electron gas. Using our own VMC energy of a uniform electron gas calculated for a system with A'^ = 284 electrons yields 
a — 338 erg/cm^, which still is much larger than the LDA surface energy. 
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where Acj) is the so-called surface dipole barrier, and Ep represents the Fermi level relative to the mean interior 
electrostatic potential. ph 
For a jellium surfaces, 

Ef ^^kl+ fi,,{n), (3.24) 

where fixci^) represents the xc potential of the uniform electron gas, and kp ~ {^i-K^nj^^^ is the Fermi momentum. 
Using Eq. ( ]3.24D one then obtains an exact general expression for the jellium work function: 



= A0--4-^,c(n). (3.25) 



Both kp and iixc{n) are known properties of the uniform electron gas and one only needs, therefore, to approximate 
the dipole barrier, which is usually evaluated in the LDA. I— , 

LDA self-consistent calculations of the work function of a jellium surface were carried out by Lang and KohnEil. 
The LDA xc potential decays exponentially outside the jellium surface and does not reproduce, therefore, the correct 
image-like asymptotic behaviour. Hence, LDA orbitals yield a too rapidly varying electron-density profile and values 
of the dipole barrier and the work function that are top high. Ncailocal calculations of the work function have been 
performed with the use of wave-function-based FHNCo and DMCe3 methods, and they indicate that the actual work 
function of jellium surfaces lies about 0.3 — 0.5 eV lower than their LDA counterparts. Nevertheless, an accurate 
determination of DMC work functions still requires longer runs and thicker simulation slabs, in order to eliminate the 
statistical fluctuations in the electron-density profiles and to avoid the presence of considerable finite-size effects. 

In the case of real surfaces the work function is obtained from Eq. ( 3.23| ), and one needs to determine both the 



dipole barrier and the Fermi level Ep. Now the bulk properties are not those of a uniform electron gas and Ep 
also needs to be approximated. In actual first-principles calculations of the work function both Acj) and Ep are 
usually obtained within the LDA. As the LDA uncertainties introduced in the evaluation of Ep tend to cancel out 
the corresponding uncertainties introduced in the evaluation of the dipole barrier, the LDA is expected to perform 
better in the evaluation of the work function of real solids than in the case of jellium. Indeed, LDA work functions of 
a variety of real solids have been shown to be in good agreement with experiment. This trend is particularly striking 
in the case of 4d metals, where the relative difference between first-principles LDA and experimental work functions 
is found to be within 5%. This is shown in Table |. 

IV. QUASIPARTICLE CALCULATIONS 
A. The surface barrier 
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The form of the surface barrier outside a solid is of great importanpSj for the interpretation of a variety of sucfai 
sensitive experiments, such as low-energy electron diffracticm (LEED)e3, scanning tunnelling microscopy (STM)l£j'l 
and inverse and two-photon photoemission spectroscopyEjO. This issue was addressed in Refs. |53| , ^ iPfiliF^I through 
the exact xc potential of DFT, i.e., the xc potential felt by Kohn-Sham electrons. A way of calculating the actual 
surface barrier felt by quasiparticle states of the many-body system is provided by the Green-function theory and, in 
particular, the-|nonlocal state-dpeendent electron self-energy 'S{r,r' ; E). 

Deisz et alS3 and White et alSB simulated the physics of the self-energy near the surface in terms of an effective state- 
dependent local potential, whose real part was interpreted as the actual potential felt by a given quasiparticle state. 



Comparing the nonlocal Hamiltonian entering Eq. (2.13) with an effective local Hamiltonian, this state-dependent 
effective local xc potential was defined by 

Viociv) ^s{r) = J dr' S(r, r', Es) *.(r'). (4.1) 



Calculations of the local potential of Eq. (LI) were reported in Ref. ^ for the lowest-lying image resonance in 
a jellium surface with = 2.07. This effective potential was found not to depend sensitively on the quasiparticle 
state and to nearly coincide with the xc potential of DFT obtained in Ref. IsJ by solving an exact integral equation, 



Eq. (3.12). Nevertheless, the calculations reported in Ref. ^ for quasiparticle states at the (111) surface of real 
Al, within 1.5 eV of the vacuum energy, indicated that the relative contributions of exchange and correlation to the 
image potential felt by unoccupied quasiparticle states are significantly different from the corresponding contributions 
to the exact Vxc of DFT. These authors showed that in the case of quasiparticle states above the Fermi level the 
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exchange part of the effective xc potential of Eq. (4.1) displays exponential decay rather than power-law behaviour, 
and concluded that even for a semi-infinite metal inclusion of correlation is essential for the correct description of the 
image potential felt by unoccupied states. 



B. Image-potential states 



An electron outside a metal surface can be trapped by the image-potential induced surface barrier, if its energy lies 
below the vacuum level and is reflected from the metal due to the lack of allowed bulk states along certain symmetry 
directions. This situation occurs when for a given wave vector parallel to the surface k|| there is a gap in the projected 
bulk band structure. Noble and transition metals usually present a surface band gap near the vacuum level Ey, 
and the long-range character of the image potential then gives rise to a Rydberg-like series of b ound image-potential 
induced states whose binding energies converge for k|| = towards Ey as dictated by Eq. ( 3.19 ). In the absence of a 
band gap, these image-states often persist as image rjesonances. 

More than two decades ago, Echenique and Pendrjfl carried out a theoretical investigation of the integrity of image 
states at metal surfaces by viewing these surface states as standing waves of an electron bouncing back and forth 
between the semi-infinite solid and the surface barrier. Image states at metal sttriaces were first detected experimen- 
tally byJLEED fine-structure analysis and then by inverse photoemission (IP)E3'L3Eil and two-photon photoemission 
(2PPE)EJo. Although the Kohn-Sham cigenfunctions and eigenvalues have no clear physical meaning, model xc 
potentials that go beyond the LDA have been used to reproduce, within IjWXt, the experimentally observed binding 
energies and effective masses of image states on a variety of metal surfacea£jE3. 

The finit|e-.lifetime of image states at metal surfaces, which is known to be mainly due to electronic excitations 
in the solidEa, is inherently outside the realm of DFT. In the framework of Green- function theory, one identifies the 
inverse quasiparticle lifetime as follows 



-2Im£;, 



(4.2) 



where Eg is the quasiparticle energy. Assuming that the cigenfunctions 4>s (r) and eigenvalues Ss of a suitable hermitian 
single-particle Hamilto nian Hp approach the actual quas iparticle w ave functions and energies, Es can be approximately 
determined from Eq. ( ^.16 ). Introduction of Eq. ( ^.16 ) into Eq. (O) then yields, on the energy-shell {Zs = 1), 



T„ = — 



2 y dr y dr' cj)lir) Iml](r, r'; e,) 0,(r). 



(4.3) 



If one further assumes translational invariance in the plane of the surface, single-particle wave functions are of the 
form 



(A,(r) = ^</>,(z)e^'^ii 



(4.4) 



and the inverse quasiparticle lifetime is obtained as follows 



dz 



dz' 4>l {z) ImE(2;, z'; k|| , Es) (t>s{z), 



(4.5) 



where 4>s{z) and now are image-state wave functions and energies describing motion normal to the surface, and A 
represents the normalisation area. 

The first quantitative evaluation of the lifetime of image states, as obtained from Eq. (|4.5|), was reported by 



Echenique et aLE23. In this calculation, the image-state wave functions (psiz) were approximated by hydrogenic- 



like wave functions with no penetration into the solid and a simplified free-electron-gas (PEG) model was used to 
approximate the electroiL^elf-energy. In subsequent calculations the penetration of the image-state wave function into 
the crystal was allowedllH^, thereby accounting for this new decay channel. The role that the narrow crystal-induced 
Shockley surface state on the (111) surfaces of Cu and Ni plays in the decay of the n — 1 image state on these surfaces 



was investigated by Gao and LundqvistliHi. In this work, the image-state wave functions were also approximated by 
hydrogenic-like wave functions with no penetration into the solid, a simplified parameterised form was used for the 
Shockley surface-state wave function, and screening effects were neglected altogether. 

The first self-coiasiatent many-body calculations of image-state lifetimes on noble and siiaplejnctals were reported 
only very recentlyll23, and good agreement with experimen tally determined decay timestSil t23 was found. In Ref. 
|l03| , the decay rate of image states was computed from Eq. (O) with the electron self-energy evaluated in the G^W^ 
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approximation. The image-state wave function (f>s{z) and energy Ss and all the eigenfunctions and eigenvalues entering 
the non-interacting Green function G'^{z, z'; k|| , e) were obtained by solving a one-particle Schrodinger equation with a 
realistic one-dimensional model potential, and the potential variation in the plane parallel to the surface was considered 
thorugh the introduction of the effective mass. 

Self-consistent calculations (ff^he key role that the Shockley surface state plays in the decajz_c)f image states on Cu 
surfaces were later carried outE£3, 
the so-called GWT approximatio 



the inclusion of short-range xc effects was investigatedt£3 in the framework of 



given by Eq. (^.19 ), but with an effective screened interaction 



In this approximation, the electron self-energy is of the GW form, i.e., it is 



W{r,r';E)=v{r-r')+ J dn J dra [w(r - n) + /..(n, ra; £;)] x(i-i, r2, ^) v(r2 - r'). 



(4.6) 



t he d ensity-response function now being the time-ordered c oun terpart of the retarded density-response function of Eq. 
(3.8). The xc kernel /ac(ri, r2; E) entering Eqs. (3^) and (4^) account for the reduction in the e-e interaction due to 



the existence of short-range xc effects associated to the image- state electron and to screening electrons, respectively. 
In Ref. 108 the exact xc kernel was replaced by that of Eq. ( 3.10| ), and it was concluded that G^VF^ calculations 
produce decay rates that are within 5% of the more realistic G"WT calculations. Although the presence of short- 
range exchange and correlation between screening electrons significantly enhances the decay probability of image 
states, this enhancement was found to be more than compensated by the large reduction in the decay rate produced 
by the presence of a xc hole around the image-state electron itself. 

Unlike image states, which are unoccupied states with energies close to the vacuum level, crystal-induced Shockley 
surface states are known to exist— neai^he Fermi level in the F-L projected bulk band gap of the (111) surfaces of 
the noble metals Cu, Ag, and Au 



111 



113 



Hence, these surface states form a quasi two-dimensional (2D) electron gas, 
which overlaps in energy and space with the three-dimensional (3D) substrate, and represent a promising playground 
for lifetime investigations. 



Photohole lifetimes of Shockley surface states were investigate i^-iji a variety of metal surfaces with the use of high- 
resolution angle-resolved photoemission (ARP) spectroscopy 



le 



Recently, the STM was used to determine 



the 



m 

118 



lifetime of excited holes at the edge of the partially occupied surface-state band on the (111) face of noble metaL 
and also to measuce.the lifetime of hot sur face -state and surface- resonance electrons, as a function of energy 
G°W° calculationaiij, as obtained from Eq. (4^), demonstrated that the decay of surface-state holes is dominated by 
2D electron-electron interactions screened by the underlying 3D electron system, and showed an excellent agreement 
with the experiment. These theoretical investigations were then extended to the case of hot surface-state and surface- 
resonance electrons in Cu(lll)EH2l, showing that, contrary to the case of surface-state holes, major contributions to the 
e-e interaction of surface-state electrons above the Fermi level come from the underlying bulk electrons and thereby 
giving an interpretation to the measurements reported in Ref. 118 



V. PERSPECTIVES FOR FUTURE RESEARCH 



Approximate density-funjCtipnal treatments of many-body effects are being applied to increasingly more complex 



problems in surface scienceE£ll. The question of their accuracy in describing these effects, e.g., in chemical reactions 
at surfaces, and how to improve them systematically is still outstanding. In the case of a jellium surface many-body 
effects are not masked by the effect of the crystal structure, so this system provides an stringent test of the accuracy 
of approximate density functionals in strongly inhomogeneous systems. In principle, approximate density functionals, 
as well as other many-body approaches, could be benchmarked for the jellium surface against diffusion Monte Carlo 
results. However, the presence of rather large statistical fluctuations (in the electron-density profile) and finite-size 
errors (in the surface energy) makes it difficult to give a clear verdict on the performance of LDA, GGA, and RPA 
on the basis of the current DMC data. There is a clear need for new DMC calculations of the work function and 
the surface energy of jellium surfaces which are performed with the use of larger system sizes, longer runs, and an 
improved treatment of finite-size effects. 

Although computationally much more expensive than DFT calculations, QMC algorithms have the important 
computational feature of being inherently parallel. The use of massively parallel computers has made the application 
of QMC methods to the investigation of real solids possible. Furthermore, QMC studies of real surfaces are now 
within reach, and a full-many body description of chemical reactions at surfaces appears to be an exciting perspective 
further down the line. 

Recently, the Green-function theory, within the GW approximation, has been shown to provide total energies of 
an infinitely extended free-electron system that are_as accurate as those obtained from DMC calculations, as long as 



the GW approximation is treated self-consistently££j. Within this approach, the Green function entering the GW 
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self-energy of Eq. (2.19) is evaluated self-consi stent ly from the Dyson equation [Eq. ( ^.12| )], the screened interaction 
being obtained as in the RPA through Eqs. (2_^)-(2^)-but with the non-interacting Green function G^{r,r';E) 



replaced by its self-consistent interacting counterpart! 



The total energy is then obtained from the so-called 



Galitskii-Migdal formulaEj. The application of this way of obtaining the total energy to the case of a jellium surface 
seems to be promising in the effort to understand the current discrepancies between DFT and wave-function-based 
estimates of the surface energy. 

Excited-state properties are still the realm of Green function methods (in the GW approximation) . In the case of 
image-potential induced states these methods have the great advantage of providing, within the same computational 
framework, a consistent description of both binding energies and lifetimes, as well as the surface barrier. So far, 
full GW calculations have been limited to the study of jellium surfaces and the (111) surface of Al, both of which 
support only image resonances. However, the binding energies and intrinsic lifetimes of these resonances are not well- 
defined and are difficult to compare with the experiment. On the other hand, currently available GW calculations 
of well-defined image states on noble metal surfaces have been carried out with the use of a model one-dimensional 
potential, and first-principles many-body descriptions that are based on a complete three-dimensional treatment of 
the band structure of the solid would be desirable. Since binding energies of image states are very sensitive to the 
surface barrier, comparison between these calculations and the available experimental data should also help to resolve 
remaining questions, such as the precise many-body origin of the image potential. 
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Surface 

Nb(llO) 

Mo(llO) 

Pd(lll) 

AG(lll) 

Ag(lOO) 

Ag(llO) 



LDA° 
4.66 
4.94 
5.53 
4.67 
4.43 
4.23 



EXP" 
4.87 
4.95 
5.6 
4.74 
4.64 
4.52 



TABLIfiJ. LDA and experimental values of the work function for low-index crystal faces of 4d metals, taken from "Ref. 
and ''Ref.H, respectively 
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FIG. 1. Contour plots of exchange-correlation hole for jellium surface for Vg = 2.07 calculated within the DFT-based RPA. 
Snapshots show deformation of the hole around the electron as the electron moves from a point well inside metal to the vacuum. 
The vertical line represents the edge of the positive background. One of the circles represent the position of the electron. The 
other circle represent the centre of gravity of the hole. Metal occupies z > 0. 
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